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ABSTRACT 



A deep ocean capsule is proposed capable o.f sensing 
pressure, temperature, flow and direction for a study of 
tides and internal waves of the oceans. One of the major 
requirements for this capsule is the detection and recording 
of unusual events. To meet this need a very low power, 
highly reliable digital detector has been designed. COS/MOS 
logic with maximiam use of MSI and LSI was used in the de- 
sign. The detector was simulated using proposed word 
lengths and fixed point arithmetic which would be utilized 
in the ensuing hardware. The fixed point simulation was 
necessary to verify the problems in truncation and round-off 
errors and to provide necessary design criteria to ensure 

I 

proper filter functioning and integral detection. 
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I . INTRODUCTION 



An unusual event detector is to be a part of a one-year 
deep ocean experiment [ Ref. 1 ] . The sensors and record- 
ing system are all digital and have low power and high 
reilJ.ability as overriding design constraints. 

The detector is to have several distinct subsystems as 
.=:J;inwn in Figure 1: a spurious data rejection gating system; 
a semi-diurnal and diurnal tide rejection filter; and a 
square and integrate fading memory detector. Upon determi- 
nation of an unusual event characterized by significant 
energy in the spectrum from 0.25 CPH to 30 CPH [Ref. 1], 
various permanent memory and accoustic telemetry links may 
•±»e initiated. The sensor signals are to be sampled at one- 
Tninute intervals and stored in a two-hour buffer memory. If 
no event is occurring, the signals are lost at the end of 
.two hours. If an event is occurring, the buffer memory is 
tJed Into the permanent memory. 

The spurious data rejection system is accomplished by 
tapping tiie buffer memory at the half-hour points and taking 
first di.fferences of adjacent signals. If the differences 
exceed a certain tolerance level, a spurious data point or 
glitch i<^ ass\amed, and the first difference is rejected and 
is replaced with an appropriate value. These five first dif- 
ference .arc then combined with appropriate weighting to 
achieve >’ u anti-tide filtering. The resulting signal is 
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Unusual Event Detector - Pressure Input is Monitored for an Unusual 

Event Which Will Trip Control Circuitry to 





then squared and fed into an integrator with a fading 
memory time constant of one-half hour. That is, the detec- 
tor will be predominantly sensitive to events that have 
accumulated over the past half-hour. 

Because the system is to operate unattended at the 
3,000-fathom level for one year, power and reliability are 
tantamount design considerations. For this reason, COS/MOS 
IC's are to be used with maximum utilization of existing LSI 
and MSI components. Based on these considerations, the 
following hardware design is proposed as a preliminary 
prototype for system evaluation. 

A logic design was formulated, and its veracity was 
analysed via a simulation using actual sensor word lengths 
and appropriate fixed point arithmetic. The simulation as 
a design tool verified an optimum design with proper 
trade-offs between system performance and power and size 
projections . 
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II. DESIGN CONSIDERATIONS 



A. SYSTEM EQUATIONS 

Although the four input variables of temperature, 
pressure, and the two current components are sampled, only 
the pressure data is sampled for unusual events in the 
tsunami spectrum of 1 to 30 CPH [Ref. 1] . A more sophisti- 
cated device could use any number of detectors on any of the 
inputs using the same temporary storage buffer memory. 

1. First Difference and Limiter 



drift, the first stage in the detector is the first dif- 
ferencer. It computes the difference between successive 
samples - Pt-1* Iri order to protect the detector from 
any spurious data which would give a false indication, a 
limiter which limits the first difference according to the 
following relation is the next step in the detector: 



2 . Tide Filter 

The tide filter is a fourth order filter with zeroes 
chosen to remove the diurnal and semidiurnal tides. Two 
sections of the form 



In order to remove the mean pressure and long-term 




«t + ^l^t-A ’^l'^t-2A 



Js (t-T)W3^(x) 



dx where 



Wj^(x) = 1 + aj^6(A) + bi6(2A) 
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with a]^, and ^21 ^2 chosen to give the required zeroes 
are combined by convolution to give the fourth order 
difference equation: 

W(z) = 1 + a]^a 2 <S (A) + jaia2 b]^b^6(2A) + 

[^1^2 a2b^6(3A) + b2^b2<5 (4A) 

With A = 1/2 hour, the five inputs to the filter will be 
the first difference 6P^ at the times 0, -30, -60, -90, and 
-120 minutes. The constants are calculated to be 
ao = 1 

ai + a2 = -3.920392 
a^^a2 t ^ 1^2 ~ ^ " ^41800 
a^b2 + a2b]^ = -3.920392 
bib2 = 1 

3 . Integration 

To obtain the power spectral density of the tide 
filter the signal is next squared and averaged over a fading 
function with a 60 minute time constant. 

^(t) = aF2'(t-A) + P2(t) a < .991666 

The only function remaining is to determine if the power has 
exceeded a preset threshold indicating significant energy in 
the tsunami band to signal an "event." 

4 . Control Functions 

The system proposal calls for a dump of the previous 
two hours of one-minute samples whenever a tsunami is 
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detected. This will ensure that an adequate history 
leading up to an event is stored. Also included in the 
logic is the circuitry to continue storing the one-minute 
samples until the second consecutive two-hour interval which 
is below threshold. Extra logic gates are provided, however, 
for a variety of controls lending flexibility to the 
detector as a "black-box" component. 

B. SYSTEM ARCHITECTURE 

The logic design for a special purpose digital machine 
begins with the system equations, block diagrams and a gen- 
eral understanding of overall system requirements. Special 
purpose structures can usefully be described by the general 
system pictured in Figure 2 illustrating the two major parts 
of a special purpose machine, the data processing and 
timing and mode circuitry. 




Figure 2. Special Purpose Structure 
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The data processing circuitry includes all input and 
output interfacing circuits, the data operators and system 
blocks, any storage of memory devices, and, in general, the 
necessary hardware to take system inputs and derive system 
outputs. The timing and mode circuitry serves much the 
same purpose for the special purpose structure as the soft- 
ware program does for the general purpose computer. The 
timing circuitry sequences operations in a set pattern from 
some initial time. The mode circuitry controls operations 
which are functions of the data being handled as well as 
time . 

In considering the best approach to the design of the 
digital package, the need for low power has repeatedly been 
stressed. For this reason, serial processing circuitry 
which has the advantage of economy of hardware at the ex- 
pense of speed is used predominantly with only slight 
modifications to meet other of the specified system 
requirements . 

1 . Memory 

The major component and power requirement of the 
detector will be the temporary data storage section or buf- 
fer memory. Investigation into this problem centered around 
core memories, shift registers, thin films and delay lines. 
Core memories have the definite advantages of nonvolatility, 
simplicity, zero quiescent power, reliability and random 
access capability. The same is basically true of the new 
thin film adaptations of» ferrite storage. Shift registers 
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are popular as short terra buffer storage, are serial in 
oper.ation, and require no external read write hardware. 

Since random access, nonvolatility of information, 
and access time are of little importance here, the most 
favorable devices are long shift registers which are compa- 
tible with the logic family. A discrete state shift register 
of length 10,000 bits at first seems to be rather prohibitive, 
but recent trends in MSI are making these long shift regis- 
ters well-suited for small size memories. Using current low 
power MOS technology as many as 64 static bits are available 
on a single chip. Since such long shift registers used as a 
recirculation memory in a serial machine need simply a 
single port input, a single port output and a clock, the 
total amount of hardware involved for the complete memory is 
very small. 

Also considered were dynamic shift registers, which 
have more than three times the number of bits per chip as 
static registers [Ref. 2] . Although the quiescent power is 
lower, the total power is greater than for static devices 
because of the transient power of continuous clocking. The 
resulting long static discrete state shift register chosen, 
therefore, was dictated by available hardware at low power, 
serial operation for simplicity, recirculation at an extreme- 
ly low rate for reduced hardware count and data accessibility, 
and low clock rate for low duty-cycle. The organization of 
data in the shift register memory is shown in Figure 3. 



15 



2 . 



Data Organization 



The type of data organization and corresponding 
arithmetic code to be used for the detector results from 
the form of the input data. All five input channels (pres- 
sure, temperature, current components and time) will be 
sixteen bit binary magnitudes. For this reason, the Data 
can be simply word organized within the samples which are 
sequentially organized as shown in Figure 3. The arithmetic 
code chosen was initially a one's complement code for simpli- 
city. Further study, however, indicated the use of a two's 
complement code was more compatible with serial processing 
and would ultimately require fewer operations and hardware. 
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Figure 3. Data Organization in Memory 
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Two additional considerations which will be clarified in 



the later section regarding the use of RCA's COS/MOS logic 
family are the use of positive AND/OR logic for initial logic 
design and the use of asychronus timing circuits. The short- 
cuts involving the use and marriage of NAND and NOR logic to 
reduce the package count is a matter of circuit juggling to 
make optimum use of the available gating packages. Asychro- 
nus or pseudosynchronus timing of all operations again 
results from the tradeoff of speed for simplicity. 

3 . Flexibility 

Two of the specifications for the unusual event 
detector are compatibility and flexibility. The input sec- 
tion of the device is therefore designed to accept a variety 

3 

of input combinations. The digital sensors for the deep sea 
in most cases are frequency modulated signals corresponding 
to the particular input variable. They are measured by 
allowing the output frequency to count a binary counter for 
some predetermined period [Ref. 3] . 

The input section is designed to accept five 
variable frequency or variable period signals. The exact 
period of measurement (from one second to four minutes) or 
the choice between frequency or time measurement is easily 
changed, using removable solder jumpers to program the input 
in a variety of fashions. 

The storage of the 120 minute samples in a static 
register provides a flexibility of output. A simple 
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clocking gate or even an external clock can be used to read 
out this inf ormation ' in almost any format. 

C. IMPLEMENTATION WITH COS/MOS 

Recent trends in semiconductor technology toward medium 
and large scale integration have made the MOS transistor 
increasingly attractive due to its ease of manufacture and 
versatility of operation. Single chip building blocks of 
rather large arrays, registers, and memories are being incor- 
porated into many logic families; of primary interest here, 
however, is the use of the MOS device in a comp.lementary 
symmetry arrangement which has the distinct advantage of 
extremely low power [Ref. 4] . 

At present the leading supplier of a complete logic line 

3 

using complementary symmetry MOS transistors is Radio Corpo- 
ration of America with its COS/MOS line. The COS/MOS CD4000A 
series of integrated circuits features a variety of options 
including special low voltage types, high reliability cir- 
cuits, various packaging styles, and the usual military or 
industrial temperature ranges. 

1 . Complementary Symmetry Circuitry 

The breakthrough that made complementary symmetry 
logic practical in medium scale integration was the ability 
to fabricate opposite polarity, but similar devices in large 
numbers on the same chip. The basic gate circuit from which 
all other logic functions are derived is shown in Figure 4. 
The usual load resistor is replaced with another FET of 
opposite polarity; thus, there is no quiescent current 
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except for the leakage currents on the order of nanoamps. 
The- pair of MOS transistors are simply switches which tie 
the output to either the supply voltage or to ground. In 
true complementary logic, the circuit for the. P devices is 
the dual network function of the N units. 




Figure 4 , Three Input Complementary Symmetry MOS NOR GATE 

The advantages derived by the complementary 
arrangement are more than just the microwatt standby power: 
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a. Simple one supply direct-coupled logic using 
only MOS transistors 

b. Large fanout capability because of capacitive 
inputs 

c. Good noise immunity due to large threshold 
voltages 

d. High speed operation because output node 
capacitance is always charging or discharging 
through the "on" unit 

2. Shift Registers Using COS / MOS 

The P and N channel units can form a shift register 
using true complementary logic with two flip-flops per bit 
for two phase operation; however, a significant reduction in 
components can be obtained by the use of (1) the excellent 
capacitance storage feature of the MOS gates, and (2) the 
two halves of the MOS flip-flop, one of which will provide 
the function of delay [Ref. 4]. As shown in Figure 5, 
opening the switches S]^ and S 2 which form the cross-coupling 
connections results in two separate storage elements, each 
capable of storing for an extended but finite period of 
time. In this state, new information can be input to the 
left gate, while the level of the right half is transferred 
to the next stage on the right. Closing the switch S 2 and 
then closing Sj^ locks the one-zero information indefinitely. 
The two phase clocking requires a rather quick transfer of 
data from the left to the right halves of a stage, but an 
indefinite delay is allowed between stages. 
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Figure 5. Basic Shift Register Stage 



3 . Power Dissipation 

Considering the total power dissipation as the svm 
of both the quiescent and transient powers, further power 
reduction can be obtained by proper design of the switching 
characteristics. Minimizing the logic level swings, using 
true complementary loads, and minimizing leakage currents as 
well as all output and load capacitances form the general 
approach to an optimum system. The total power P dissipated 
in a complementary set-reset flip-flop is given by 



P = + P3 

where Cq is the node capacitance, Vq 
f is the operating frequency, and Pg 



is the supply voltage, 
is the standby power. 
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For low switching power, the supply voltage Vq must be 
reduced to the lowest possible value. At a Vq of -4 volts, 
and an output capacitance Cq of 23 pf the frequency depen- 
dence is shown in Figure 6. Note, also, the large increase ' 
in quiescent power with temperature [Ref. 4] . 



MICROWATTS 
FOR A 
SINGLE 
FLIP-FLOP 
STAGE 




0 5 10 15 20 kHz 



Figure 6. Power Dissipation as a Function of Repetition Rate 

Since only a fixed number of operations are completed 
each minute in the detector, the clock rate will have no ef- 
fect. The power economy is gained by reducing the total num- 
ber of redundant operations. The repetition rate can go 
down to DC; however, the longest allowable high state of the 
clock input of the 64 bit shift registers can be no more than 
100 microseconds. For this reason, the clock frequency was 
chosen to be 19.360 kHz which provides a good margin of 
safety . 
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III. HARDWARE IMPLEMENTATION 



The logic design formulated in this section is organized 
according to the logic blocks diagrammed in Figure 1 of the 
Introduction. All logic functions are described and defined 
in the text and the figures. The later sections appraise 
the cost, size and power requirements of the resulting logic 
packages . 

A. CLOCK AND TIMING 

The first section of the detector, which sequences all 
operations and controls the functions of the device, is the 
clock and timing circuitry. Table I, which may be consider- 
ed as the program, lists all sequences and defines the 
timing signals originating in the COS/MOS counter circuitry. 

Since the COS/MOS inverters have nanowatt standby power 
and microwatt operational power, an excellent low power clock 
can be designed using the forward gain of an inverter with a 
suitable feedback arrangement. A good amplifier section for 
the oscillator is formed with an inverter stage and feedback 
resistor R-|^ large enough (>10M) so that the attenuation and 
phase of the feedback are not appreciably affected. The 
resistor biases the output at about half the supply voltage. 
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Figure 7. Crystal Oscillator Circuit Using 
RCA CD4007 Inverter 

The crystal controlled phase shifting network is added 
in Figure 7 with provision for trimming to the design fre- 
quency. Using this circuit, stability should be on the 
order of 5 ppm and total supply current from a supply voltage 
of 4 volts is 50 microamps [Ref. 6 ] . 

Frequency division of the crystal oscillator is 
necessary due to duty cycle considerations which place the 
logic clocking frequency at 19,360 Hz as shown in Figure 8. 
One CD4004 provides four flip-flops in an asychronus counter 
to divide the 311.6 kHz crystal frequency down to the clock- 
ing frequency of 19,360 Hz. The bit counter consists of one- 
half of a CD4013 dual flip-flop plus a CD4022 eight bit 
divider/decoder to aid in decoding the timing signals. The 
word counter is another CD4022 connected to divide by five 
with a CD4001 gating package to reset the counter. 
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The sample counter is a CD4004 seven bit counter 
connected to count to 121, and ten other logic packages 
provide the necessary reset and decode functions. The 
second counter is a similar arrangement connected to reset 
at 60 using the same type logic packages. The first stage 
in the second counter is the recirculation counter; it 
divides by two to give a frequency of 1 Hz. 
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Figure 8. Clock - Bit and Word Counters Plus Decoders to 

Derive Timing Signals 
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TO MINUTE COUNTER 
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TABLE I. 



TIMING 



These Boolean expressions define sequentially the timing 
3j_gn3ls referred to on the later diagrams. The gating cir 
0-jjitry for these expressions is shown in Figures 10 /II /and 12. 



Sec(0)«R(0) RECIRCULATE BUFFER MEMORY 

Sec(0)»R(0) FIRST DIFFERENCE AND LIMITER 

S(L)»E ]3 = Sec(0)*R(0) [s(0)+S(30)+S(60)+S(90)+S(120)3 
S(H).Ej 3 = Sec(0)*R(0) [S(1)+S(31)+S(61)+S(91)+S(121)J 



Ed-W(O) 

S(H).Ed[w( 0) + W(2j] 
S(H)-Ed-W(1) 

S (H)*Ej)- W(2) 

S (H)*Ed* W(2)« B (15) 
Sec(l)* R(0) 

Ed*W(2) [b( 8) to B(15j]* 



^sum* W (0) 



Clock 16 bit delay register 
and adder 

Shift Limiter register 

Check difference and Set 
latch if good. Set sign 
change latch if necessary 

Shift present difference 
to temporary storage. 

Change sign. 

TEMPORARY STORAGE 

Shift in 8 bit difference 

Shift out difference 

TIDE FILTER 



Sec (2)*R(0) 

Esum = (2) •R(0) 

Esum[w(0) B(5) toB (15)|*+W(1) Multiply by -3.920392 

Es^[w(0) B(7)toB(15^*+W(l) Multiply by 5.841800 

E_^* S (1) + S(121)] • Add differences 1 and 121 

|?(1)*B(4) to B(12]]* 

Es^ [s (1)+S (12irj»W(l) ‘B (12) Set sign latch 
^s;am + ♦ W (1) • B ( 0^ Shift accumulator 
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TABLE I. Continued 



Sec(3) R(0) 

Egg = Sec(3)*R(0)* [s(0) 
Egq'Wd) • B(0) 

Esq«W(0) 

Sec (4) • R(0) 

Eint = Sec(4)«R(0)*S(0) 

to Bds]] 

Eint-Wd) 



SQUARE 
to S(75j]* 

Shift accxjmulator in sum 
and check bit for "one" 

Shift input store and SQUARE 
accumulator 

INTEGRATE 

Shift in the squared value 
Shift registers 



Sec(59)»R(0) INPUT SECTION 



Ein = Sec(59)« 


R(0)« S (121) 






Ein«W(0) 


* 


Shift 


in pressure value 


Ein'Wd) 




Shift 

reset 


in temperature, 
pressure counter 


Ein«W(2) 




Shift 

reset 


in current^, 
temperature 


Ein«W(3) 




Shift 

reset 


in current^, 
current^ 


Ein*W(4) 




Shift in time, reset 
current^ 


* B(i) to B(j) 


= B(i) + B(i+1) + 


. . . + B ( j ) 



The decoding circuitry for these signals are shown 
together in Figures 10,11/ 12. Usually they will be 

located on the same board as the device they control. 
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Sec (0) 
R(0) 



RECIRCULATE MEMORY 




* 




FIRST DIFFERENCE 



Clock 16 bit delay 
register and adder 



Shift Limiter register 

Check difference and 
reset latch if good. Set 
sign change latch if 
necessary 

Shift present difference 
to temporary storage 

Change sign 



S?H) 

W(2) 

B(H) 




TEMPORARY STORAGE 
Shift 8 bit difference 



Sec (2) 

R(0) 

W(0) 




Shift out difference 



Figure ‘ 10 . Timing Decoding Gates 
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TIDE FILTER 



Sec(2) 

R(0) 



W(l) 



W(0) 

B(H) 



B(5) 

B(6) 

B(7) 



W(l) 



W(0) 




Multiply by -3.920392 



Multiply by 5.841800 



B(7) 

B(H) 



S(l) 

S(121) 

B(H) 

W(l) 

^sum 

B(4) 

B(L) 

B(4) 

B(5) 

B{6) 

B(7) 

B(H) 

B(0) 

B(l) 

B(2) 

B(3) 




Figure 11. Timing Decoding Gates 
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SQUARE 



Sec (3) 
R(0) 

[S(0)toS(7) ] 



W(l) *8(0) 



W(0) 



W(l) 

Sec (4) 

R(0) 

S(0) 

W(0) 

B(H) 

B(0) 

B(l) 

B(2) 




Shift accumulator in sum 
and check bit for "one" 



Shift input store and SQUARE 
accxjmulator 

Shift Registers 



Shift in the sguared value 



INPUT SECTION 



Sec (59) 

R(0) 

S(121) 

W(0) 



W(l) 

W(2) 

W(3) 



W(4) 




Shift in pressure value 

Shift in temperature, reset 
pressure counter 

Shift in current^j, reset 
temperature 

Shift in current^, reset 
current^ 

Shift in time, reset current^ 



Figure 12. Timing Decoding Gates 
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B. PROCESSING CIRCUITS 



1 . Input Section 

The five inputs consisting of variable frequency 
signals are used to either count or gate five sixteen bit 
binary counters. The fifth counter is provided with the 
additional circuitry so that it will count time in minutes. 

On command from the system clock each of the five counters 
sequentially shifts its value into the system memory and, 
except for the time channel, resets for the next counting 
interval. Figure 13 shows only the pressure and time channel 
because all five inputs are the same except for the gating 
on the minute counter. Since the two-hour sequencing signal 
for the detector is obtained from the time counter, it must 
reset at 65,520, a multiple of 120 minutes, so that the 
two-hour sequence is maintained. 

All of the counters are made from CD4004's. Three 
counters for each input provide the necessary sixteen bits. 
Two eight bit parallel to serial shift registers convert 
the counter value into serial form for the system's buffer 
memory. Fourteen additional packages provide the gating and 
mode control circuitry. 
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Figure 13. Input Circuitry - Each Output Can be Wired to 

Either Gate or Count the 16 Bit Counters. Thus, 
the Mean Frequency or Mean Period of the Input Is 
Shifted to the 2 Hour Buffer Memory Once Each 
Minute 
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Figure 13. Input Circuitry - Each Output Can be Wired to 

Either Gate or Count the 16 Bit Counters. Thus, 
the Mean Frequency or Mean Period of the Input Is 
Shifted to the 2 Hour Buffer Memory Once Each 
Minute 
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2 . System Memory 



The system memory consists of a discrete state static 
shift register of 9680 bits in a recirculation loop storing 
the past 121 minute samples organized sequentially by words 
and bits. The bit, word and sample counters provide for 
addressing the data; further division of the clock frequency 
into seconds and minutes provides all necessary timing sig- 
nals for the rest of the detector. The timing signals are ob- 
tained by coincidence gating across the stages of the asychro- 
nus binary counters. In order to reduce the transient power, 
the data is recirculated through the register length only once 
per minute. On this one pass, clocking and gating read out 
the ten pressure words needed for computation of the five 
first differences. The 9680 shifts averaged over a minute's 
time results in an average clocking frequency of 150 Hz. 

From data similar to the graph of Figure 6, the transient 
power used in the memory should be about 100 microwatts; 
however, the 64 bit shift registers are still in the devel- 
opmental stage and this power estimate is meant only as a 
comparative estimate. 

In order to obtain the 9680 bit length in COS/MOS 
logic, 152 logic chips must be series connected. The exact 
9680 bits required is made of 151 sixty-four bit TA5989 
registers and one CD4006 sixteen bit register. The total 
memory package should dissipate about 0.3 mW and fill a 
vol\ime of 0.34 cu.ft. on six standard 2.5 by 5.0 inch logic 
cards . 
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3. First Difference and Limiter 



First differences are calculated once a minute on 
the five sample pairs S(0), S(l); S(30), S (31) ; . . ,S (120) , 
S(121) by delaying the first sample in a register, comple- 
menting the next and adding the two serially. Since the ori- 
ginal input number is allowed to overflow the input counters , 
the signed result must be corrected if the sign is inconsis- 
tent. Gating also checks the difference for unreasonable 
values before allowing them to pass on to temporary storage. 

The sixteen bit delay and all the eight bit storage 
requirements in the first differencer can be obtained from 
various connections of the CD4006, an eighteen stage regis- 
ter. Two CD4015A dual four stage parallel output registers 
provide the necessary connections used in the limiter to 
sense and correct the difference value. One-third of a 
TA5963 triple serial adder is used for the adder and a 
CD4013 dual flip-flop performs the complement function. 

Twelve additional packages as listed in Table II provide 
all the gating and control logic. 

4 . Temporary Storage 

Since only correct first differences are allowed 
through the limiter, the temporary storage must be capable of 
nondestructive read-out when the tide filter calls for inputs. 
This is a very simple matter of circularly shifting the data 
to the right. Only the 8 most insignificant bits need to be 
stored, because the average first difference is on the order 
of 4 bits. The validity of this assumption was verified by 
simulation. 
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Figure 15. First Difference and Limiter 




Figure 16. Temporary Storage - Five 8 Bit Registers Store 

the Previous Differences in Case the Present One 
Is Rejected by the Limiter 
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5. Tide Filter 



The five differences temporarily stored are simply 
multiplied by appropriate weighting values and summed to ob- 
tain the filtered output. Since there are only three dif- 
ferent constants, the differences are gated to the proper 
section through serial switches and the sum of these values 
is stored in a register which also serves as the input 
register to the squarer. 

Since one of the weighting constants is unity, only 
two multipliers are required. They simply shift and add the 
input the appropriate number of times and proper clocking 
keeps the result in the same sixteen bit format. 

The constant -3.920392 can be approximated with the 
number -3.920410 or 2^ + (2~4 + 2"6 +2"10 + 2"H) . This 
constant, when multiplied by the input number of sixteen bits 
gives a possible maximum length result of 29 bits. Through 
simulation and device consideration, however, the input 
difference will rarely exceed eight bits and the approxima- 
tion to reality of the constants makes the significance of 
the lower five bits questionable. Truncation to sixteen 
bits before addition, then, will have no effect. Similarly, 
5.84180 is approximated by 5.842797 or 2^ - (2^ + 2~3 + 2 “^ 

+ 2-9) . 

The maximum of the sum of the five weighted differ- 
ences will be approximately the sum of the magnitudes (1+4 
+6+4+1) times the average input value. Limiting the 
maximum first difference to seven bits or 128^0 requires 
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eleven whole number bits in the accumulator leaving four 
fractional bits plus the sign bit. To multiply by the con- 
stant one, the input need only be left-shifted four bits to 
account for the change in the binary point. Clocking and 
gating accomplish the necessary four bit delay. 

The analytic two's complement code used throughout 
the detector requires no special handling in the multipliers 
however, since the eight bit differences are to be extended 
to sixteen bit precision words in the accumulator, care must 
be taken to pack the additional eight created bits with 
zeros or ones depending on the sign of the original number. 
The sign latch in Figure 18 tests the sign from the 1st and 
121st difference and ensures that the input is held at the 
proper level if there are no incoming bits . 

All of the timing signals referenced on the figures 
are defined in Table I. The Figures 10, 11 and 12 diagram 
how the timing signals are derived through gating and 
decoding the frequency divider outputs. 

The total number of twenty- five delays in the two 
multipliers is provided by four CD4015 dual four bit regis- 
ters. Three more TA5963 triple serial adders perform the 
summing functions, with the result stored in an additional 
two CD4015's. The two complementors are similar to the one 
used in the first difference and diagrammed in Appendix A. 
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W(1)+W(0) • [B{5)toB(15) ] 
R(0) •Sec(2) 




TO CLOCK ALL 
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W{3) 



j. 

RESET 
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. AND 
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MULTIPLIER X -3.920392 APPROXIMATED IN BINARY TO 26 BIT 

PRECISION BEFORE TRUNCATION TO 16 
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Figure 17. Multiplier X 5.841800 APPROXIMATED IN BINARY 
TO 26 BIT PRECISION BEFORE TRUNCATION TO 16 
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18. Accumulator. Each of the Weighted Differences Are Added into the 

Accumulator maintaining 4 Fractional Bits 
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6 . Squaring 



The squaring package observes the value stored in 
the sum accumulator of the filter and shifts and adds it into 
the output storage. The square circuitry would handle nega- 
tive as well as positive numbers if care were taken to pro- 
vide all the necessary ones in the extra negative bits. A 
simpler method to implement, however, is to simply rectify 
the input. By testing the sign bit in the sum accumulator 
before each operation, negative numbers can be complemented 
bit by bit. This method, as shown in Figure 19., gives an 
error of one in the least significant bit because it takes 
the one's instead of two's complement. The error is 
insignificant . 

For each bit the input is shifted and if the current 
bit is a 1, the shifted value is added to the result. The 
sum in the output of the filter contains a sign and four 
fractional bits; in order to make the detector adaptable to a 
variety of dynamic ranges, the tap on this register is made 
adjustable. The input to the squarer can only be eight bits 
in length, but the adjustable tap gives a range of the most 
sensitive eight bits from 1/32 to 15 to the maximum values 
between 16 and 65,536. Simulation later showed the most 
probable tap to be at the binary point. 
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Figure 19. Squares 8 Bit Value from Tide Filter by Shifting One Bit, Testing for 
SQUARER ”1" and Adding Incrementally to the Accumulator 




7. Integrator 



Since the squared output is to be integrated over a 
decaying time constant of sixty minutes, the maximum avail- 
able storage must be about thirty times the maximum value to 
ensure against overflow. To stay within the sixteen bit pre- 
cision, the input squared value must be truncated to the 
eleven most significant bits. The decay factor of .991666 
is accomplished by delaying seven bits and then subtracting. 
Taps are accessible on the delay register to vary the time 
constant of integration. Coincidence gating on. the storage 
register synchronized to the clock determines the value in 
the tsunami detector and provides the control logic with its 
go/no-go signal. 

The necessary storage for the squarer and integrator 
can all be formed from the CD4006 eighteen stage register. 

A CD4015 dual four bit register provides for the multiplica- 
tion in the integrator along with two of the summers left 
over from the triple serial adder used in the first differ- 
encer. The complementer again comprises a dual flip-flop 
CD4013 and nine additional gating packages complete the 
squarer and integrator. 
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8 . Storage Control 

This section of the detector is designed with flexi-^ 
bility as the prime consideration since the actual device used 
for long term storage will be an independent system. Whether 
a tape system or core storage is used the only output required 
of the detector is a yes or no answer every two hours quali- 
fied by considerations of its previous answers. A critical, or 
event, is defined whenever the value in the Integrator regis- 
ter exceeds the threshold set by the sense tap. A store 
command is given under the following conditions : 

a. A critical is detected and the past n consecutive 
hours have not been critical; 

b. A critical is not detected and the past two hours 
were critical. 




Figure 21. Control - Three Flip-Flops and a Counter 
Provide Logic to Determine if an Unusual 
Event Should Be Signaled 
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C. PACKAGE CONSIDERATIONS 



•The COS/MOS line offers sufficient complexity per chip 
to keep the overall package count low, and also, enough var- 
iety to provide for optimum utilization of the available 
circuitry. The total package requirements for each of the 
sections of the detector are listed in Table II, with a 
cost estimation in Table III. 

1 . Space and Cost 

The size and cost predictions of the detector are 
based on an experimental set-up which will use dual in-line 
packaging on a general purpose bread board arrangement. 
Standard Logic Inc. offers a line of logic cards which have 
various combinations of fourteen and sixteen pin DIP sockets 
and all mounting and interconnecting hardware. Computer aided 
wiring and interconnections complete the system. Using the 
Standard Logic boards from 20 to 30 DIP packages depending 
on the number of sixteen and fourteen pin packages, can be 
arranged on each 2.5 by 5.0 inch card. Figure 22 shows the 
logic arrangement and overall dimensions of this experimental 
set-up. 

For reliability the detector should be made with 
flat-packs on printed circuit boards with all solder con- 
nections. The space requirements will be approximately 40 
to 60 percent of the experimental arrangement due to the 
elimination of the DIP sockets and the interconnecting 
wiring . 
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TABLE II. 



PACKAGE REQUIREMENTS 



FUNCTION 

CLOCK 



INPUT 



MEMORY 



FIGURE 


NUMBER 


LOGIC < 


CIRCUIT 


8 


2 


CD4022 


Divide by 8 Counter/Decoder 




1 


CD4004 


Seven Stage Counter 




1 


CD4013 


Dual D Flip-flop 




1 


CD4009 


Hex Inverter 


9 


11 


CD4023 


Triple 3-Input NAND Gate 




4 


CD4001 


Quad 2-Input NOR Gate 




2 


CD4009 


Hex Inverters 




2 


CD4004 


Seven Stage Ripple Counters 




1 


CD4013 


Dual D Flip-flop 




25 


total 




13 


15 


CD4004 


Seven Stage Ripple Counter 




10 


CD4014 


8-stage Parallel in/Serial out 




5 


CD4011 


Quad 2-Input NAND Gate 




4 


TA5872 


Dual J K Flip-flop 




2 


CD4023 


Triple 3-Input NAND Gate 




1 


CD4000 


Dual 3-Input NOR plus Inverter 




1 


CD4013 


Dual D Flip-flop 




1 


CD4009 


Hex Inverter 




39 


total 




14 


151 


TA5989 


64-Stage Static Shift Register 




2 


CD4011 


Quad 2-Input NAND Gate 




1 


CD4006 


18 stage Static Shift Register 




154 


total 
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TABLE II. Continued 



FUNCTION FIGURE 

FIRST 15 

DIFFERENCE 



TEMPORARY 16 
STORE 



MULTIPLIERS 17 



NUMBER LOGIC CIRCUIT 



5 CD4011 

2 CD4013 

2 CD4015 

2 CD4016 

2 CD4012 

2 CD4009 

1 CD4006 

_1 ^TA5963 

17 total 



Quad 2-Input NAND Gate 
Dual D Flip-flop 
Dual 4-Stage Parallel out 
Quad Bilateral Switch 
Dual 4-Input NAND Gate 
Hex Inverter 

18 Stage Static Shift Register 
Triple Serial Adder 



3 CD4006 18-Stage Static Register 

3 CD4011 Quad 2-Input NAND Gate 

3 CD4016 Quad Bilateral Switch 

2 CD4009 Hex Inverter 

_2 CD4023 Triple 3-Input NAND Gate 

13 total 



4 CD4015 

3 TA5963 

2 CD4011 

2 CD4013 

_1 CD4023 

13 total 



Dual 4-Stage Parallel out 
Triple Serial Adders 
Quad 2 -Input NAND Gate 
Dual D Flip-flop 
Triple 3-Input NAND Gate 
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TABLE II. Continued 

FUNCTION FIGURE NUMBER LOGIC CIRCUIT 



SUM 18 

ACCUMULATOR 



SQUARE 19 



INTEGRATE 20 



2 CD4011 Quad 2-Input NAND Gate 

2 CD4015 Dual 4-Stage Parallel out 

1 TA5963 Triple Serial Adder 

1 CD4023 Triple 3-Input NAND Gate 

1 CD4000 Dual 3-Input NOR plus Inverter 

1 ^CD4002 Dual 4-Input NOR Gate 

8 total 

2 CD4006 18-Stage Static Register 

2 CD4011 Quad 2-Input NAND Gate 

1 CD4001 Quad 2-Input NOR Gate 

1 CD4009 Hex Inverter 

1 CD4016 Quad Bilateral Switch 

1 TA5963 Triple Serial Adder 

8 total 

3 CD4015 Dual 4-Stage Parallel out 

2 CD4011 Quad 2-Input NAND Gate 

1 CD4013 Dual D Flip-flop 

1 ^TA59 6 3 Triple Serial Adder 

7 total 
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TABLE II . Continued 



FUNCTION 


FIGURE 


NUMBER 


LOGIC CIRCUIT 




STORAGE 


21 


2 


CD4011 


Quad 2-Input NAND 


Gst0 


CONTROL 




1 


CD4004 


Seven Stage Ripple 


Counter 






1 


CD4013 


Dual Flip-flop 








1 


CD4009 


Hex Inverter 








5 


total 






TIMING 


10 


6 


CD4011 


Quad 2-Input NAND 


Gute 


GATES 


11 












& 12 


2 


CD4001 


Quad 2 -Input NOR Gate 






2 


CD4023 


Triple 3-Input NAND Gate 






2 


CD4025 


Triple 3-Input NOR Gate 






1 


CD4002 


Dual 4-Input NOR Gate 






1 


CD4012 


Dual 4-Input NAND 


Gate 






2 


CD4009 


Hex Inverter 





16 total 
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Figure 22. Physical Configuration 



I 







TABLE III 



COST EVALUATION 



Low Voltage Series 


CD4000A 3 to 


15 volts 




Flat Pack 


-55°C to +125° 


C 






Number 


Cost Each . 


Total Cost 


CD4000AK 


2 


$ 7.00 


$ 14.00 


CD4001AK 


6 


7.40 


44.40 


CD4004AK 


19 


14.75 


280.25 


CD4006AK 


7 


18.00 


126.00 


CD4009AK 


12 


11.00 


132.00 


CD4011AK 


31 


7.40 


222.00 


CD4012AK 


3 


7.65 


22.95 


CD4013AK 


9 


10.05 


90.45 


CD4014AK 


10 


20.15 


201.50 


CD4015AK 


11 


20.25 


222.75 


CD4016AK 


6 


10.80 


64.80 


CD4022AK 


2 


20.95 


41.80 


CD4023AK 


19 


7.65 


145.35 


CD4025AK 


2 


7.65 


15.30 


TA5872 




10.00 


30.00 


TA5963 


7 


14.00 


100.00 


TA5989 


151 


20.00 


3020.00 


TOTAL 






$4752.55 


Similar arrangement in 

Dual In Line Plastic Packages 


Approx- 


$2000.00 



-40° to +75°C 
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2 . Power 



The actual quiescent power required varies 
considerably from chip to chip and is very temperature 
dependent. In most cases, RCA specifications include a 
"typical" and "maximum" value. The estimates given assxme 
the "typical" measurement and an ambient temperature of 
25®C. The actual power at the ocean floor, therefore, 
should be less than the estimate. Also the transient power 
is of no significance in the digital processing circuitry 
because the maximum number of operations of any of the cir- 
cuits is still on the order of 160 per minute. This gives 
an average repetition rate of only 3 Hz when averaged over 
the one minute period. For this reason only the transient 
power of the clock, the decoder and the buffer memory are 
included. 

The total power drain of 0.5 mW required for the 
whole detector system is almost three orders of magnitude 
smaller than a comparable logic package of conventional cir- 
cuitry. Compared to the power requirements of the sensors, 
accoustic links, and other functions of the capsule, the 
power for the detector is negligible. The power for one 
year's operation should be obtainable from a single 2 amp- 
hour storage battery. The power used to run the crystal 
oscillator should be well regulated and be run from a sepa- 
rate mercury battery for voltage stability. All of the other 
detector circuits can be supplied from lead-acid cells due to 
the excellent voltage compatibility of the COS/MOS line. 
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If power drain were not so important, more than a 
two to one reduction in package count could be effected by 
rebuilding the memory out of conventional MOS registers . 
The power required would increase almost ten-fold. Future 
developments in the COS/MOS line toward LSI might even 
reduce the size as much with no sacrifice in power. 
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TABLE IV. 



POWER REQUIREMENTS 



LOGIC ( 


CURRENT 


STATIC POWER 


NUMBER 






Microwatts 




CD4000 


Dual 3-Input plus Inverter 


0.01 


2 


CD4001 


Quad 4-Input NOR Gate 


0.01 


6 


CD4002 


Dual 4-Input NOR Gate 


0.01 


2 


CD4004 


Seven Stage Ripple Counter 


1.5 


19 


CD4006 


18 Stage Static Register 


0.05 


7 


CD4009 


Hex Inverter 


0.05 


12 


CD4011 


Quad 2 -Input NAND Gate 


0.005 


31 


CD4012 


Dual 4-Input NAND Gate 


0.005 


3 


CD4013 


Dual D Flip-flop 


0.025 


9 


CD4014 


8 Bit Parrallel in/Serial out 


4.0 


10 


CD4015 


Dual 4-Stage Parallel out 


4.0 


11 


CD4016 


Quad Bilateral Switch 


0.1 


6 


CD4022 


Divide by 8 Counter/Decoder 


1.5 


2 


CD4023 


Triple 3-Input NAND Gate 


0.005 


19 


CD4025 


Triple 3-Input NOR Gate 


0.01 


2 


TA5872 


Dual J K Flip-flop 


0.025 


4 


TA5963 


Triple Serial Adder 


1.5 


7 


TA5989 


64-Bit Static Shift Register 


1.0 


151 


TOTAL i 


STATIC POWER 


268.8 Microwatts 


TRANSIENT POWER: CLOCK 


270 






MEMORY 


110 





TOTAL TRANSIENT POWER 380 Microwatts 
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IV. COMPUTER SIMULATION 



A preliminary step in the design of the system components 
was the evaluation of the required functions through computer 
simulation. The general algorithm has previously been ana- 
lyzed [Refs. 1,5]. System evaluation of the actual logic sys- 
tem should be as close to actual hardware circuitry as 
possible with fixed point numbers and identical computer logic 
simulation . 

A. COMCOR Cl 5000 

For actual checkout of some of the functions in hardware, 
the digital patching section of the Cl 5000 which is located 
in the Electrical Engineering Computer Lab, was used. The 
complementer for the two ' s complement arithmetic was the 
first function wired. The second logic block that was test- 
ed extensively in hardware before design was finalized was 
the multiplying scheme. Since a single multiplier used all 
32 of the computer's general purpose flip-flops and a major 
portion of the NAND gates, this method was applicable to only 
the system blocks; another computer had to be chosen for the 
overall system evaluation. 

B. CDC 160 

With the consideration that the simulation should retain 
as much simularity to the actual hardware as possible, the 
CDC 160 computer, also located in the Electrical 
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Engineering Computer Lab, was chosen. The distinct 
disadvantage, later to become obvious, was the twelve bit 
word limitation in trying to check out a sixteen bit system. 
In most areas of the program and logic flow, there were no 
differences between the designed sixteen bit system and the 
twelve bit simulation; however, the checkout of the multi- 
pliers required further investigation. At this point the 
168 arithmetic unit was used as peripheral equipment to 
allow for double precision. 

1 . Adaptation to Fixed Point 

A magnetic tape was available from the I.G.P.P. 

containing two records of pressure variations at the 2000 
fathom level in an area near San Clemente. The data had 
been taken at four minute intervals over a four day period 
with both a Vibratron and a Hewlet-Packard sensor. The data 
was linearly packed to simulate the required one minute sam- 
pling rate, and this prepared tape was used [Ref. 5] . A 
major portion of the program listed in the Appendix was re- 
quired to read the BCD tape and convert it into pure binary 
form. Several problems due to word length limitations were 
solved by reading in only the four least significant decimal 
digits and converting them to binary through a table look-up 
scheme. Simple testing of the resultant value for consis- 
tency removed any ambiguity in operations with the numbers 
except when a change in the fifth BCD digit occurred as noted 
on Figures 24 and 27. This method was incorporated into the 
actual design also because the input counters were allowed 
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to overflow any number of times. Figures 23 and 26 show a 
record of the tape data compared to the truncated version 
used on the 160, 

2 , Precision Considerations 

In order to avoid going to double precision and using 
the external arithmetic device, which would no longer bear 
any^ relation to the hardware being simulated, careful study 
of the data revealed several compromises. In most cases a few 
of the most significant bits could be dropped because the 
particular data did not reach the maximum excursion designed 
for in the sixteen bit system. This, combined with truncat- 
ing one or two of the less significant bits, reduced the word 
length to the required twelve bit format. No significant 
degradation in overall performance was noted from these com- 
promises except in the case of the multipliers. Here, twelve 
bit precision allows for an error of .002 which is not insig- 
nificant since the least-count resolution of the input sen- 
sors is on the order of .0005. The multipliers were tested 
for proper functioning in twelve point individually, but, to 
allow for system evaluation and to obtain usable results, 
double precision was used. 

The only other notable difference between the 
simulation and the design was the use of one's as opposed to 
two's complement arithmetic. Although one's complement re- 
quires the additional operation of end-around carries to 
correct the sum on subtraction, this is handled in the ma- 
chine; therefore, on the software level, compensations were 
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made so that the simulation was the same as the two's 



complement on the hardware level. 

3 . Simulation Results 

The program operated on the input data samples in the 
same format as the proposed system, sequentially taking one 
sample value at a time and carrying it through to increment 
the final integrator storage. Because of the limitation of 
only 4K of storage in the 160, the input data had to be buf- 
fered into the routine 120 values at a time and the output 
unloaded at 1500 values per run. This incremental method of 
printing output on the Cal-Comp plotter adds an extra line 
at each beginning point. 

Since the Hewlet-Packard sensor had a very smooth 
output compared to the Vibratron, it was used for all pre- 
liminary runs to test the tide cancellation of the filter 
section. Once smooth output was obtained from the filter by 
adjusting the precision and values of the weighting constants, 
the square and integrate routines were perfected on the 
Vibratron data. Whether the perturbations detected on the 
Vibratron tape were noise or actual events is of no concern 
to the detector's operation but is a problem of sensor 
design . 

Tracing the three interesting "events" detected on 
the Vibratron tape back to the original data, indicates 
pressure perturbations on the order of 4mm. Considering the 
tidal variations as the major noise source of the input data, 
and these 4mm. variations as signal, the overall improvement 
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in signal-to-noise ratio is on the order of 2. 3*10^ or 
67 db. A more detailed software analysis of system perfor- 
mance, detection thresholds, power spectral densities, and 
evaluation is included in the Refs. 1 and 5. 
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PRESSURE FLUCTUATIONS 
FROM VIBRATRON DATA 
TAPE IN BCD FORMAT 
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Figure 23. Vibratron Sensor Output 
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PRESSURE VARIATIONS 
IN BCD FORM FROM 
DATA TAPE 
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Figure 26. Hewlet Packard Sensor Output 
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TIME (MINUTES) 
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Figure 28. Hewlet Packard Sensor Data 
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Figure 29. Difference of Precision in Multipliers 
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V. CONCLUSIONS 



The unusual event detector implemented with COS/MOS 
logic should provide a very useful and needed package to the 
growing number of tools used in oceanographic research. The 
detector can be used as an independent component in a larger 
system, or, due to its system structure design it can be 
used as the basis for a rather complex package. Several sub- 
systems for ascent and descent control, telemetry, data 
storage, and sensor stabilization are actually needed to com- 
plete an underwater instrument capsule. The clock and timing 
circuits of the detector could be easily extended to cover 
other operations and simple gating of the memory package and 
addition of more memory cells, could make it a general buffer 
memory for the other functions attached to the system. 
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APPENDIX A. COS/MOS 
LOGIC PACKAGES 



Lead numbers are not given since package styles vary. 




CD4001 QUAD 
TWO INPUT NOR GATES 



CD4025 TRIPLE 
THREE INPUT NOR GATE 




CD4011 QUAD 
TWO INPUT NAND GATES 



CD4023 TRIPLE 
THREE INPUT NAND GATE 
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DATA 

RESET 

DATA 

RESET 




FOUR STAGE REGISTER 
WITH SERIAL IN PARALLEL OUT 



EIGHT BUFFERED INPUTS 

CD4021A 8 STAGE SYNC 
PARALLEL INPUT/SERIAL OUT 



Vdd 



Vss 



Vc Vc 




EIGHT DECODED OUTPUTS BILATERAL SWITCH 



Vq 

Vi 

^c 

Vq 

Vi 

Vss 



CD4022 DIVIDE BY EIGHT 
DIVIDER DECODER 
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t 



CLOCK 



RESET 




CD4009 HEX INVERTERS 





















7 STAGE 

RIPPLE COUNTER 









































SEVEN BUFFERED OUTPUTS 



V 



SS 



V 



CC 



CD4009A 7 STAGE 
RIPPLE BINARY COUNTER 



^L 




CD4006 EIGHTEEN STAGE 
STATIC SHIFT REGISTER 



CD4013 DUAL D 
FLIP-FLOP 
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THE COMPLEMENTOR CIRCUIT USED SEVERAL TIMES IN DETECTOR 
DESIGN CAN BE FORMED FROM ONE CD4013 DUAL FLIP-FLOP AND 
ONE CD4011 QUAD TWO INPUT NAND GATE 





THE SERIAL SWITCH CAN BE FORMED FROM ONE CD4011 QUAD 
2 INPUT NAND GATE OR ONE-HALF OF A CD4016 QUAD BILATERAL SWITCH 
PLUS AN INVERTER 
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I 



I 



COMPUTER PROGRAM FOR SIMULATION 
OF THE UNUSUAL EVENT DETECTOR 



CrEWINDJ 1000 

iJf 



MEMORY LOCATION AT BEGINNING 
OF FUNCTION BLOCK 








1520 




.700 


- 


IYIU-LiTIFIjY 




SQUARE 






i 


► 


1550 


1 


7300 


INTEG^E^ 



4740 



( STORE IN y 

V T ,..,BUFFERyp„LL ^ 



PLOT ON 
CAL COMP 






INITIAL 
OUT BUFFER 

__I 



) 



( SHIFT ) 

RF.GTq .TF.R,^ 



1600 



76 




0000 

0000 

0300 

0300 

0300 

0300 

2100 

o4oo 

7101 

1560 

2167 

0720 

6203 

7101 

1260 

okoo 

U167 

0300 

2030 

ouo 

7101 

1560 



OPERATICKAL RCUrDIES FOR VARIOUS ENTRY 
AND EXir-OUTPUT OPTICNS 



nop REMOVE ALL HALTS EXCEPT AT REWIND, OUTPUT, AND 

nop COMPUTATia? COMPLEH'E 

nop 

nop 



1 dm TO READ OUT INPUT DATA 



Jfi 



Idi 

sbn 



1567 1600 enter 
1242 0704 enter 
1245 1260 enter 
1327 3330 enter 



1150 

0705 

1250 

7250 



Idn 

sti 



7250 2100 1 dm truncated, version 



o4oo 

7101 jfi 
1560 



Idd 

ls 3 

Jfi 



TO READ CUT FIRST DIFFERENCE 



1412 0570 enter 

1425 7200 enter 



o4oo 

7240 
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0000 

0101 

o4oo 

1000 

1050 

1100 

1150 

3350 

lUUo 

1520 

1550 

l6oo 

1700 

l?4o 

2000 

4000 

7000 

7200 

7300 

7U00 



CSWERAL LESTING OF TSUNAMI DETECTOR PROGRAM 



to 


0100 


TEMPORART STORAGE LOCATICKS 


to 


0377 


CCNVERSICN miBERS+ BINART TO OCTAL 


to 


0770 


SHIFT REGISTER STORAGE 


to 


1047 


GENERAL INITIALIZE AREA 


to 


1070 


ZERO INPUT ARRAY AND BEGIN TO READ 


to 


1322 


READ 120NUMBERS lOTO 2000 to 3777 


to 


3336 


BEGIN INPUT CCfTVERSiai TO BINARY 


to 


1426 


DIFFERENCE ROUTINE 


to 


1533 


PREPARE FOR DOUBLE PRECISia? 


to 


1550 


CALL MULTIPLY SUBROUTINE 


to 


1570 


SUM THE WEIGHTED DIFFERENCES 


to 


1637 


SHIFT REGISTER IN 0400 to 0570 


to 


175^ 


SQUARE ROUTINE 


to 


1777 


INTEGRATE ROUTINE 


to 


3777 


INPUT BUFFER STORAGE 


to 


6777 


OUTPUT BUFFER STORAGE 


to 


7161 


MULTIPLY ROUTINE 


to 


7224 


CLIPPER ROUTINE 


to 


7375 


CAL-CCMP PLOTTER ROUTINE 


to 


7777 


RS-022 FLEX AND OCTAL READ ROUTINES 
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LISTING OF STORAGE LOCATICKS SET ASSIDE FROM 



oooo 



0077 



TO 



0001 into mult for 2 and 4 

0002 temp storage of shifts 2 

0003 ts of output of mult for 2 and 4 

0004 ts value of p entering mult for 2 and 4 

0005 

0006 ts of shifts 6 

0007 ts 



0010 

0011 ts for shift register routine 

0012 

0013 

0014 

0015 

0016 

0017 



0020 

0021 


0400 


location of 


0022 


0401 




0023 


0600 


location of 


0024 


0601 




0025 

0026 

0027 







shift register words one 
shift register words two 



0030 

0031 difference 1 

0032 difference 31 

0033 differmce 6l 

0034 difference 91 

0035 difference 121 

0036 

0037 



0040 

0041 mult diff 1 

0042 mult diff 31 

0043 mult diff 6l 

0044 m mult diff 91 

0045 mult diff 121 

0046 

0047 
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0050 

0051 

(X)52 

0053 

0054 

0055 

0056 

0057 

0060 

0061 

0062 

0063 

0064 

0065 

0066 

0067 

0070 

0071 

0072 

0073 

0073 

0074 

0075 

0076 

0077 



indlgit + correction 

running sxjm first word 

manning sura 2nd word 

address of indigit being operated on 

4th digit of converted nos. and how many over 4 

7000 store 

count of nonzero and not 20 digits 
address of first of 4 indigits to be used 



reserve all remaining for rs-022 
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0000 


7101 

71^05 


Jfl 


LOAD FLEX TAPE USING rso22 


1000 


7500 

ll6i 


exc 


REWIND TAPE AND BEGIN PROBLEM AT START 




2200 

4 ooo 


Idc 


initialize constants 




iaoo 

1561 

4100 

1747 


stm 

stm 




1010 


4100 

1771 

o 4 oo 

4 ooo 


stm 

Idn 

std 






4 ooi 

7700 

7101 

1050 


std 

hit 

Jfl 


zero the flex read command for safety 
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1050 2200 


Idc ZERO INPUT ARRAY AND BEGIN TO READ 


2000 




4100 


sjSm 


1056 




o4oo 


Idn 


4100 


stm 


2000 




5701 


aoB 



io6o 


3600 

3770 


sbc 




6706 


njb 




2200 

1777 


Idc 




4063 


std 




7700 


hit halt when array is zero 




7101 




1070 


1100 
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BEGIN INPUT COTVERSiai TO BINARY 



1150 2200 


Idc 


0100 




4100 


stm 


3262 




o 4 oo 


Idn 


4 o 6 i 


std 


4062 


std 


0300 





1160 


888 

000 




1163 


5463 


aod 




2163 


Idi 




6005 






o 4 oo 


Idn 




4057 


std 


1 


1170 


7101 

1220 


m 




5457 


aod 




0704 


sbn 




6003 






7101 

1163 






7700 


hit 


1 


1200 


2200 

2000 


Idc 




4100 

1122 


stm 




4063 


std 




o 4 oo 


Idn 




4057 


std 




7101 

1050 


6 t ± 



zero 6l and 62 



step 63 by one and read new digit 
if new digit is zero, step 3 and count 
zero the zero count 



test for 4 consecutive zeroes 
halt when out of input digits 

initialise constants 



go to zero array and read more indigits 
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2163 

0720 

6oo4 

5466 

7101 

1163 

2066 

0704 

6205 

0400 

4o66 

7101 

1163 

4o64 

o4oo 

4o66 

2063 

0300 

0704 

4067 

7101 

1260 



Idl bego;n to ccnvert to binary 

sbn 

zjf if 20 step 4 
aod 

Jfl step 66 1 and go read another 
Idd 

sbn if it was 20 



pjf are past 4 real nos 
Idn 

std zero 66 
Jfi 

read more 

std 

Idn 

std zero 66 



Idd 

sbn 

std store in 67 address of 1st of 4 digits to be converted 
Jfl 
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2167 

5200 

0100 

4o6o 

2061 

3160 

4o6l 

0500 

2200 

0200 

5060 

2062 

5160 

4062 

0420 

5100 

1262 

5600 

0200 

6oo4 

5467 

7101 

1260 

0300 

2061 

1200 

7000 

4o64 

2061 

5464 

4100 

o4oo 

2o64 

0110 

3062 

4100 

0600 

7700 

7101 

3350 



Idl BEGIN TABLE LOOK UP 
adc 

std 

Idd 

adl 

std add new digit to srunning stun 



Idc 

rad add 200 to indigit + corr. nos 

Idd 

adi 

std 

Idn 

ram 



sbc 

zjf test to see if we have read 4 nos 
aod 

Jfi go back to I260 and convert more 



Idd separate nos into 2 words 
Idc 

std 

Idd 

sbd store 12 least significant bits in shift 
stm register 0400 to 0570 



Idd 

ls3 

add 

stm store most significant bits in 0600 to 0770 
hit halt input to binary is complete, now filter 
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2200 

0400 

4020 

2200 

o6oo 

4022 

0430 

4024 

2020 

o6oi 

4021 

2022 

o6oi 

4023 

2120 

3521 

4124 

6203 

7101 

1372 

2122 

3523 

6012 

6305 

o4oo 

4124 

7101 

1410 

2124 

1600 

7000 

4124 

7101 

1410 

2122 

3523 

6013 

6206 

o4oo 

0300 



Idc DIFFERENCE ROUTINE 

std Initialize constants in 20,22,24 

Idc 



std 

Idn 

std 



Idd begin each difference here 

adn 

std 

Idd 

adn 

std 

Idi difference word 1, least significant digits 

sbi 

sti 

Pjf 

Idi for positive difference of word one 
sbi 

njf 



Idn for pos word 1 and pos word 2, zero the difference 

sti 

Jfi 

Idi for pos word 1 and neg word 2, comleroent the sign 
see 

sti 

jn 



Idi 

sbi 

POf 

Idn 



for negative difference of word one 



for neg word 1 and neg word 2, zero the difference 
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4124 

flOl 

I4l0 

2124 

1600 

7000 

4124 

0300 

2020 

3600 

0570 

6010 

0436 

5020 

0436 

5022 

0402 

5024 

7101 

1340 

7101 

7200 



fftl 

Jfi 

Idl 

see 

stl 



for neg word 1 andpos word 2 , eomplement the si 0 i 



Idi test for eompletion of all 5 differenees 
she 

zjf 

Idn step the dlfferenee ecunter by 36 

rad 

Idn 

rad 

Idn 

rad 

jfl go back to 1340 and calculate another difference 
Jfi go to spike removal routine 
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m 






1 




.0450 

4026 

0601 

4027 

2126 

6206 

2200 

7777 

4127 

7101 

1455 

o4oo 

4127 

2026 

0740 

6206 

2026 

0602 

4026 

7101 

1442 

2200 

7000 

4oo4 

2200 

4505 

4o42 

2200 

7770 

4o45 

2200 

2570 

4o44 

2200 

0015 

4o45 

2200 

0002 

4025 

2200 

0000 

4o24 

7101 

1520 



Idn PREPAKE FOR DOUBLE PRECISICN 

Std 

adn 

std 

Idi for negative nos. store ones in 12 most significant 

zjf bits 

Idc 

sti 

Jfi 

Idn for positive nos, zero roost significant bits 

sti 

Idd 

sbn 

Idd 

adn 

std 

Idc 

std store multiplication routine access location. 



Idc load constant -5.919921875 

std 

Idc 

std 

Idc load constant t5» 84179700 

std 

Idc 

std 

Idc load constant +1,000000000 

std 

Idc 

std 

Jfl 

jump to call ARITH routine 
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1520 0101 pta CALL MULTIPLE SUBROUTINE 

7004 Jpi 

0051 multiply difference 1 

0050 
0024 
0046 

0051 multiply difference 31 

0032 
1530 0042 

0050 

0031 multiply difference 6 l 

0034 

0044 

0052 

0031 multiply difference 91 

0036 
1540 0042 

0054 

0031 multiply difference 121 

oo4o 

0024 

0056 

7101 

1550 Jvmg? to sximm the weighted differences 
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1 






^7 

3051 

3053 

3055 

3057 

0500 

0300 

0500 

Uioo 

Uooo 

5701 

3600 

6770 

6203 

7101 

1600 

7700 

7101 

7300 



Idd SUM THE WEIGHTED DIFFERENCES 

add 

add 

add 

add 



stm 

store the result in UOOO 

aob 

she 



Pjf 

jfi jvm^) to the shift register 

hit halt when UOOO to 6770 Is full of output 

Jfi 
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2200 

0570 

4 on 

^11 

o6oi 

4012 

2111 

4112 

2011 

0701 

4 on 

5600 

o 4 oo 

6612 

0300 

0500 

2200 

0770 

4013 

2013 

0601 

4 oi 4 

2113 

4 n 4 

2015 

0701 

4013 

3600 

0600 

6612 

7101 

1150 



Idc SHIFT REGISTER 0400 tO 0571 

std 

Idd 

adn 

std 

Idi shift data to right one sample in 0400 to 0570 
sti 



Idd 

sbn 

std 

she 

PJ1> 



Idc shin other half in 0600 to 0770 

std 

Idd 

adn 

std 

Idi 

sti 



Idd 

sbn 

std 

sbe 

pjb loop arcimd again 

go back to convert another input nos to binafcy 
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4050 

7101 

1700 

2200 

7000 

4oo4 

2200 

0100 

4o54 

0500 

0300 

o4oo 

4051 

4053 

4055 

4057 

2450 

6302 

4050 

0101 

7004 

0031 

0050 

0050 

0052 

0031 

0052 

0054 

0056 

2057 

7101 

1560 



std ENTER THE SQUARE ROUTINE 

Jfi 



Idc square routine 
std 

Idc scale factor (o4oo divides ty 2) 

std 



Idn 

std 

std 

std 

std 

led 

njf 



zero most significant digits 



rectify 



pta enter ARITH subroutine 
Jpi 

multiply squares input 



multiply scales output 



Idd load into 5^ the squared and shifted results 

Jfi 
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0500 

2100 

6770 

4o4o 

2057 

1200 

7770 

OIU 

0110 

4o57 

2o4o 

0111 

0277 

4o42 

2o4o 

5442 

6202 

o4oo 

3057 

7101 

1770 

4100 

4000 

5701 

0701 

4100 

1742 

7101 

1565 

0111 

0110 

0105 

1200 

3777 

4on 

0500 

1771 

1771 

2200 

6770 

4100 

1742 

0500 



nop INTEGRATICN ROUTINE 
1dm 

std input first value 

Idd 

Ipc 

ls6 

ls5 

std divide squared input by 8 

Idd 

ls6 

Ipn 

std 

Idd 

sbd 

Idn takes care of negative zero 

add 

Ofl 



stm store the integrator output in the out. buffer memory 

aob increment memoiy one 

sbn 

stm 

jump back into program 



ls6 corrections to CPICCMP PLOTTER ROUTINE 

ls5 for scaling 

ls2 

Ipc 

std 



Idc corerections to initial section 
stm 
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1 

I 




A 







2200 

ivooo 

UlOO 

7311 

7500 

44oi 

oUoo 

4oio 

2100 

4ooo 

4ou 

7101 

7320 

2011 

5410 

4012 

6205 

7101 

7340 

7101 

7330 

7402 

2012 

0701 

4012 

66o4 

7401 

7101 

7350 

7401 

5412 

6702 

7101 

7350 



Idc CAL-CCM> PLCdTEH ROUTIKE 



stm 

exc 

Idn 



Idxa 

std 



Idd compute difference from last point 

sbd 

std 

Pjf 



write plus X one 

Idd 

sbn 

std 

pjb 

-X 

Jfi 



write minus x one 

aod 

njb 
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7350 o 3 oo 
7koU 

5500 


write minus y 

aom step input location one 


7311 

36 cx) 


sbc 


6770 

6314 

7700 


njf 

hit halt, output buffer is empjiy 


7360 2200 
4 ooo 


Idc initialize storage number 4000 


4100 


stm 


1561 

4100 


stm 


1747 

4100 

1771 


stm 


7370 7101 
1600 




2011 

4010 

7101 

7310 


Idd 

std 
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0750 

0720 

0670 

o6Uo 

o6io 

0560 

0530 

0500 

0450 

0000 

0144 

0310 

0454 

0620 

0764 

0130 

0274 

o44o 

o6o4 

0000 

0012 

0024 

0036 

0050 

0062 

0074 

0106 

0120 

0132 

0000 

0001 

0002 

0003 

0004 

0005 

0006 

0007 

0010 

0011 

0000 



conversion for 1st 3 octal digits of 4th Indiglt 



conversion for 3rd Indlgit 



conversion for 2nd indlgit 



conversion for 1st indlgit 
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0001 

0005 

0005 

0007 

0011 

oou 

0015 

0017 

0021 

0000 

0000 

0000 

0000 

0000 

0000 

0001 

0001 

0001 

0001 

0000 



conversion forlast 2 octal digits of 4th indlgit 



conversion for last 2 octal digits of 3rd indigit 
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! 

I 



3 



1 



AKITH 

closed subroutine for control of the CDC l 68 Arithmetic Unit. 



PUBPOSE: 

This generator routine provides the capability to implement 22 -bit 
fixed point arithmetic sequences using a three address code. Both 
fractional and integer arithmetic are provided. An indirect addressing 
option is included. 



USAffiJ: 

LOAD the ARITH subroutine anywhere in memory. (Space required l^lg) 
Store the load address in any convenient low-core directory cell (OOXX); 
CALLING SEQUENCE IS; 



• 

PTA 

JPI XX 
OP 
A 
B 
C 
OP 
A 
B 
C 



OP Codes 

01 - add 

02 - subtract 

50 - multiply fractional 
31 - multiply integer 

4 0 - divide fractional 

4 1 - divide integer 

A end B: Operands 

C : Result 



JNXT) 



NOTES; 



1 , Any nuniber of operations may be called for in sequence 
prior to return to normal 160 coding, 

2 . Operands must be in standard format for 22 -bit fixed- 
point arithmetic (l^s complement with allowance for 
sign bit in each word). Least significant half-word 
in EVETI numbered cell. Most significant half-word in 
ODD numbered mate, 

5. FOR normal direct addressing operands must be referenced 
only by even numbered address. 

4 , Use of odd addresses will provide indirect addressing. 

5, Cells 70-77 are used as tengjoraries by the program. 

6 , Timing - Approximately 900 microsecs, per operation. 
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rrifite 




0100 

0101 

0102 

0105 

0104 

0105 

oio6 

0107 

0110 

0111 

0112 

0113 

on4 

0115 

on6 

0117 

0120 

0121 

0122 

0123 

0124 

0125 

0126 

0127 

0130 

0131 

0152 

013>3 

0134 

013>5 

0156 

0137 



ARITH 



SUBROOTINE FOR DOUBIE PRECISICK 



OPERATICNS USING THE CDC I68 







rem 




sr-ARITH 






3rem 




control routine for I68 xmit 


0602 








mlc - 8/30/64 




adn 


02 


4 o 77 




std 


77 


form address of 1st argument 


0470 


inlt 


Idn 


70 


set to fill OP-A-B-C in 


4076 




std 


76 


70-71-72-75 


0504 




Icn 


04 




4075 




std 


75 




2177 


loop 


Idi 


77 


transfer arguments 


4176 




stl 


76 


modify as required 










if odd( indirect) 


0201 




Ipn 


01 




6011 




zjf 


incr 




2075 




Idd 


75 




o6o4 




adn 


04 




6006 




zdf 


incr 




2176 




Idi 


76 




4202 




stf 


02 




2100 




Idi 







0000 








4176 


sti 


76 




5477 Incr 


aod 


77 




5476 


aod 


76 




5475 


aod 


75 




6517 


nzb 


loop 




2077 


Idd 


77 




4007 


std 


07 




2200 scode 


Idf 




set to generate sel code 


5323 




3523 


3323 for plus 


4270 


stf 


sel 


3365 for minus 


2070 


Idd 


70 


3303 for molt 


0701 


sbn 


01 


3302 for divide 


6015 




setout 




0701 


sbn 


01 




6011 


2jf 


sub 
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oi4o 

Ol4l 

6 i 42 

0143 

0144 

0145 

0146 

0147 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

0160 

0161 

0162 

0163 

0164 

0165 

0166 

0167 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 



0730 

6304 

0521 dvd 

5257 

6106 

0520 nul 

5254 

6103 



o44o sub 

5251 

2071 setout 

4262 

0602 

4247 

2072 

4257 



0602 

4245 

2070 

0740 

6335 

6007 

0471 dvl 
4247 



0602 

4235 

6312 

6471 hop 
0471 dvf 
4240 
0602 
4225 



2071 
4235 
0602 

4223 

2072 dlv 
4232 
0602 

4224 



sbn 


30 


njf 


mul 


Icn 


21 


raf 


sel 


nzf 


setout 


Icn 


20 


raf 


sel 


nzf 


setout 



Idn 


4o 


raf 


sel 


Idd 


71 


stf 


fl 


adn 


02 


stf 


tl 


Idd 


72 


stf 


f2 



adn 


02 


stf 


t2 


Idd 


70 


sbn 


4o 


nJf 


act 


zjf 


dvf 


Idn 


71 


stf 


f2 



adn 


02 


stf 


t2 


nzf 


div 


zjb 


init 


Idn 


71 


stf 


fl 


adn 


02 


stf 


tl 



Idd 


71 


stf 


f2 


adn 


02 


stf 


t2 


Idd 


72 


stf 


f3 


adn 


02 


stf 


t3 



set output addresses A-B 



these addresses OK for 
plus, minus, multiply 



process divide integer 



HOP links REPEAT to INIT 
process divide fractional 



common for divide 



100 



'm 



0210 

.0211 

0212 

0215 

0214 

0215 

0216 

0217 

0220 

0221 

0222 

0225 

0224 

0225 

0226 

0227 

0250 

0251 

0252 

0255 

0254 

0255 

0256 

0257 

0240 

0241 

0242 

0245 

0244 

0245 

0246 

0247 

0250 

0251 

0252 

0255 

0254 

0255 

0256 

0257 

0260 

0261 



2171 




Idi 


71 


sense sign of dividend 


6205 




Pjf 


pz 




0500 


nz 


Icn 






4071 




std 


71 




4072 




std 


72 


fill neg zero for DD 


6lo4 




nzf 


act 




o4oo 


pz 


Idn 






4071 




std 


71 




4o72 




std 


72 


fill pos zero for DD 


7500 


SlC*t/ 


exf 






0000 


sel 






activate 


7312 




out 


fl 




0000 


tl 








7311 




out 


f2 




0000 


t2 








2070 




Idd 


70 




074o 




sbn 


4o 




6307 




njf 


setln 


skip unless divide 


7305 




out 


f3 




0000 


t5 








6l04 




nzf 


04 




0000 


fl 








0000 


f2 








0000 


f3 








2073 


setln 


Idd 


73 


set input address C 


4220 




stf 


f45 




0602 




adn 


02 




4203 




stf 


t4 




4207 




stf 


t5 




7214 




inp 


f45 




0000 


t4 








2070 




Idd 


70 




0730 




sbn 


30 




6103 




nzf 


05 


skip unless MUF 


7207 




inp 


f45 




0000 


t5 








2107 




Idl 


07 




1200 




Ipf 






7700 






7700 


sense valid op code 


^64 


repeat 


zjb 


hop 


no valid - repeat 


7007 




Jpi 


07 


valid “ exit 


0000 


f45 








0000 




end 
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